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Abstract 

I have performed scaling tests using quenched spectroscopy of a family 
of fermion actions which have a hypercubic kinetic term, gauge connec- 
tions built of smeared links, and an anomalous magnetic moment term. 
These actions show improved rotational invariance compared to the stan- 
dard Wilson action and to the tadpole-improved clover action. Hyperfine 
splittings are improved compared to the standard Wilson action (at the 
level of a factor of three in the lattice spacing), and are about the same 
as for the tadpole-improved clover action. 



1 Introduction 



I report on tests of a family of fermion actions for lattice gauge theory simula- 
tions, which are designed to improve scaling of hadron spectroscopy. While very 
expensive to simulate, they appear to more than repay their computational cost 
with better scaling behavior. 

The actions which I tested were inspired by the fixed point (FP) action 
program |j for fermions[|[ ||, ||, ||, 0, |[ g, [HJ. However, they are not FP 
actions. They are Wilson-fermion like, in the sense that they have four com- 
ponent spinors on all sites and there are no manifest symmetries which protect 
the bare quark mass from being additively renormalized. 

The new features which I tested include: 

1) a hypercubic kinetic energy term. Each fermion in the action communi- 
cates with 3 4 — 1 = 80 nearest neighbors. This term improves the hadron dis- 
persion relation compared to that from actions using a standard on-axis nearest 
neighbor coupling. 

2) Gauge connections built of very fat gauge links, links built by averaging 
the fundamental link variables over a local region and re-projecting them back 
onto the gauge group. This construction removes short distance fluctuations 
from the correlators during the simulation process, rather than attempting to 
divide them out at the end. This results in observed very small additive mass 
rcnormalization of bare quantities, and I conjecture that all perturbative cor- 
rections to observables are considerably reduced. To include very fat links in 
a simulation with dynamical fermions using known technology might be quite 
expensive, but in quenched simulations the cost is minimal. 

I also tested a complicated lattice anomalous magnetic momentum term 
in which the quark and antiquark are not fixed to the same site. Some kind 
of term is needed to correct the lattice free quark magnetic moment and to 
improve meson and baryon hyperfine splittings. However, my tests show that 
the standard "clover" term, suitably normalized, improves scaling as much as 
the complicated Pauli term needed to satisfy the FP equations. 

The scaling tests in this paper are a bit non-standard (when compared to 
other studies of spectroscopy) and deserve some explanation. I am interested in 
scaling tests which are uncontaminated by extrapolations in volume or to the 
chiral limit. Thus I compare only simulations in fixed physical volume. Any 
volume would do, and so I choose a small one simply because these actions 
are expensive to simulate. To set the scale I use a gluonic observable because 
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these are quenched simulations. As a choice of gluonic observables, one has T c , 
the critical temperature for deconfinement, the string tension a or the Som- 
mer radius r$. Glucball or torelon measurements arc just too costly. Of these 
observables, the ones associated with the potential (a and ro) require a fit to 
a function V(r); the choice of the fitting function can affect the results. At 
very coarse lattice spacing this problem becomes more serious. Thus I use T c 
to set the scale, and do simulations on lattices of fixed size L = 2/T c . With 
~ 440 MeV and y/a/T c ~ 1.60, the scales are T c = 275 MeV, L = 1.45 
fm, and 2tt/L = 860 MeV. Lattice spacings aT c = 1/2, 1/3, 1/4 correspond to 
a = 0.36, 0.24, and 0.18 fm. 

All the tests are performed at fixed physical quark mass (defined cither 
by interpolating lattice data to a fixed value of m^/rrip or to a fixed value of 
tn^/Tc). At very coarse lattice spacing, and with heavy (though still relativistic) 
quarks, scaling violations from conventional actions are very large. Thus only 
modest statistics are required to identify improvement-or lack thereof-compared 
to them. 

Besides spectroscopy, I measure the meson or baryon dispersion relation. In 
lattices of fixed physical volume set by T c , the physical momenta corresponding 
to the different allowed lattice modes are multiples of T c , ap — 2irri/L or p = 
irT c ri, if L = 2/T c , and one can compare data with different lattice spacings at 
the same physical momentum. Wilson and clover fermions at < 6.0 (using 
the Wilson gauge action) exhibit bad scaling or rotational invariance violations. 
The new actions are rotationally invariant even at a = 0.36 fm. 

The outline of the paper is as follows: In Section 2 I describe the new features 
of these actions. Section 3 is devoted to scaling tests, and I make some tentative 
conclusions in Section 4. I describe the (new) FP gauge action used in these 
simulations in the Appendix. 



2 Ingredients of the Actions 

2.1 Hypercubic kinetic term 

The fermionic free field action has the generic form 




(1) 
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It is constructed by finding some free FP action M, 0, [|, El. We begin with a 
continuum action for fermions which has no doublers and is chirally symmetric. 
We construct an action on a coarser distance scale A' by iterating the FP 
equation 

(AT 1 , = H>, „' 

(2) 

where f2 is a blocking kernel and k is a tunable parameter. We select a blocking 
kernel, iterate the RGT to find a fixed point action, and then tune parameters in 
fl to make the action maximally local. We have used a factor-of-two rescaling in 
which f2 is restricted to a hypercube: Ojj is nonzero only if j = i±/i, i±/x±i/, . . . 
i±p±v±\±a, and fly = c{{i — (i — (i — j)^, (i — Each site 

communicates to 3 — 1 = 80 neighbors. 

There are many good parameterizations, resulting in fairly local FP actions. 
However, one ultimately wants to use these actions in simulations, and the 
action must be somehow truncated. There are a number of (subjective) criteria 
to select a good RGT, based on the properties of the truncated action (which the 
RGT does not know about): a good dispersion relation, E(p) — \p\ out to large 
\p\ with no complex roots, good free-field thermodynamics, P = l/3erT 4 even 
at large discretization, etc. We tuned the RGT to optimize these criteria, and 
our choice is given by c(l, 0, 0, 0) = 0.03, c(l, 1, 0, 0) = 0.01, c(l, 1, 1, 0) = 0.005, 
c(l, 1, 1, 1) = 0.0025, and k = 44.0. 

The couplings of the FP action for our RGT (for massless fermions) fall off 

exponentially with r = x 1 ^. The largest entries at distance r — 2 are 

at location x = (±1,±1,±1,±1). The smallest truncation which accurately 
reproduces the main features of this FP action is to an action which sits on a 
hypercube-that is, the free field action is to be parameterized with five nonzero 
A's and four nonzero p's, corresponding to each of the nonzero offsets. An 
example of a dispersion relation for this hypercubic action is compared to the 
Wilson action in Figs, [j] and E| We show both branches of the hypercubic 
action's dispersion relation; all roots are real. The non-truncated FP action has 
a perfect dispersion relation E = \p\ for all p. 

For massive fermions, we need an action which is on a renormalized trajec- 
tory (RT) for some RGT (with the mass scaling at each step with the ratio of 
lattice spacings). To reach the RT, one can begin with an action which has 
a very small mass but is otherwise close to a FP action, perform a series of 
blockings, and follow it out. 
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Figure 1: Dispersion relation E(p) vs \p\ for mo = for the hypercubic action. 

One complication with this procedure is that an action which is local for 
small mass can block into an action for large mass which is not short range. To 
avoid this, I take an RG transformation whose parameters are functions of the 
mass and tune the parameter(s) to insure a local action at each blocking step. 
The resulting A's and p's are smooth functions of the mass. Again, the dispersion 
relation for hypercubic approximations to RT actions are well behaved out to 
large \p\. 

Plots of the variation of the parameters-or the tables of numbers correspond- 
ing to them-are by themselves not very useful for calculation. Rather than give 
them, I will immediately present a simple parameterization of the couplings. 

My choice of parameterization is to regard all the p's and A's as linear 
functions of the bare mass and to determine A(0, 0, 0, 0) for positive bare mass 
by solving the dispersion relation for it: defining 

R = 2po(l, 0, 0, 0) + 12po(l, 1, 0, 0) + 24po(l, 1, 1, 0) + 16po(l, 1, 1, 1) (3) 
A s = 2A(1, 0, 0, 0) + 12A(1, 1, 0, 0) + 24A(1, 1, 1, 0) + 16A(1, 1, 1, 1) (4) 
D = 8A(1, 0, 0, 0) + 24A(1, 1, 0, 0) + 32A(1, 1, 1, 0) + 16A(1, 1, 1, 1) (5) 
the pole of the propagator (A 2 (p) + pfi(P) — at Pn = (* m cb 0, 0, 0)) is at 

A(0, 0,0,0) = -A s cosh(m - 1) - R sinh(m ) — D. (6) 

My linear parameterization runs over < m < 0.32, which is the useful range 
of mass for light hadron spectroscopy for lattice spacings aT c — 1/2 to 1/4. 
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Figure 2: Dispersion relation E(p) vs \p\ for massless Wilson fermions. 

The big problem in parameterizing an approximate action is that the additive 
renormalization of the quark mass for g 2 ^ drives the critical quark mass m c q 
negative, outside the region where the solution of the FP equation is meaningful. 
With Wilson fermions, and one free parameter (k) this is not a problem (one 
just tunes k above k c without otherwise altering the action) but with these 
complicated actions there is not a clear cut way to proceed. In principle, all 
the parameters in the action should depend on properties of the gauge field (for 
example, on the local value of the plaquette). To circumvent this, I make an 
arbitrary choice: I assume that the parameters continue to vary linearly with 
the bare mass, and I determine A(0, 0, 0, 0) by continuing the low-mass limit of 
Eqn. |g| to negative mass. 

A(0, 0, 0, 0) = -A s - D- Rm . (7) 

Since m c q ~ —0.4 at (3 c (N t = 4) and moves toward zero at bigger /3's, A(0, 0, 0, 0) 
is basically a linear function of mo- 

This parameterization is not appropriate for studying charm with this action 
(where arriQ = 2 to 4 depending on the lattice spacing). The parameters of the 
action on the RT are smooth functions of mo, so one could probably construct 
a more complicated (polynomial or exponential) fit to them. 

Table | gives the parameters of the linear fit. 
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Table 1: 


Linear parameterization of the couplings of the hypercubic action: 


\(xq,X\, 


X2,X 3 ) = X a ~ 


f m X b , po(x Q ,x 1 ,X2,x 3 ) 


= Pa + m p b . 




offset 


Xa 


x b 


Pa 


Pb 





2.256756 


-0.9863 






10 


-0.1122 


0.0741 


-0.1464 


0.1300 


110 


-0.0323 


0.0271 


-0.0329 


0.0303 


1110 


-0.0144 


0.0141 


-0.0101 


0.0096 


1111 


-0.0072 


0.0076 


-0.0035 


0.0033 



2.2 Very Fat Links 

Measurements of pure gauge observables (the potential or glueball masses) suffer 
from noise arising from the short distance fluctuations of the gauge fields. A 
good cure for this problem has been known for many years: define new link 
variables which do not couple to the UV sector of the lattice variables, and 
which have the same IR properties as the original variables. An example of 
such a variable is an APE-blocked link (TTJ 

V? +1 {x) = (l-a)V?(x) + a /6j2(V"(x)V-(x + l ))V;(x + ^ 

+ V v n (x-i>)^(x-P)V;(x-P + il)) (8) 

(with V®(x) — Ufi(x) and V™ +1 (x) is projected back onto SU(3)). It is also 
known that for best results, both a and the maximum number of blocking steps 
N should increase as the lattice spacing decreases. 

Fermions also suffer from bad UV behavior, and their symptoms include 
the breaking of flavor symmetry (for staggered fermions), large additive renor- 
malization of the bare mass (for Wilson fermions), and large renormalizations 
of currents (for any kind of fermion). The tadpole improvement program [[l2f 
was originally designed to estimate or compute these large UV effects and sub- 
tract (or divide) them out during the conversion from the lattice calculation to 
continuum number. 

However, recent evidence suggests that it may in some cases be better to 
remove the UV fluctuations directly from the simulations. This evidence is the 
partial restoration of flavor symmetry breaking for staggered fermions by replac- 
ing the link by an V = 1 APE-blocked link, as shown by Ref. jl3| and (with 
a slightly different averaging) by Ref. |pfl . These authors restrict themselves 
to N — 1, presumably because they wish to use their actions for simulations 
with dynamical fermions. However, if one is interested in quenched simulations, 
one can APE-block to any desired level, with tiny overhead, simply by pre- 
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computing and storing the APE-blocked links. Then if N > 1 improves UV 
behavior, one is free to use it. 

It is easy to understand why the fat links suppress UV fluctuations p5[ . 
Each term in the action in coordinate space 

Cl = Ya £ ■■U li {x + y)...il){x + z) (9) 

can be expanded as a power series in g 

£i = h £ w ^ X ) T ■ ■ ■ C 1 + tsa^/ifr + y + «/2A) 

-liagfA^x + y + a/2fif ...).., ip(x + z) 



which in momentum space becomes Li — L\ + L\ with 



(10) 



4 = if / ^(P) r I ^(k + q- p)A^k)e^ +a ^ k ^(q)e t " z (11) 



$ = ^ I P , q $(p) t Lm ^ + fc2 + 1 - p)Mh)Mh) 

e i(y+a/2jX){kt+k 2 ) q prq\ e iqz _ 

(12) 

Smearing the link over a distance ro makes the replacement 
or 

M k ) -> J2 H^(k)A u (ky(y +a w-v k (u) 

fx . V 

(and a similar formula for £|) where the form factor is H^(k) = ^2 r h )lu (r)e lkr . 
Essentially any smearing function suppresses the vertex at k > 7r/ro. In the 
language of Ref . Q , tadpoles contribute beyond their naive strength because 
the UV divergence of the gluon loop compensates for the a-dependence of the 
vertex; smearing suppresses the coupling of the fermion to high momentum 
gluons. 



2.3 Nonlocal Pauli Term 



Lattice fermions have a magnetic moment which is anomalously small due to 
lattice artifacts. One can parameterize the vertex through the interaction of a 
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fcrmion with an infinitesimal magnetic field B: the pole in the propagator will 
be at E — mo + B/2rriB, where tub is the so-called magnetic mass. We write 
the momentum-space interaction term as 

^(p)»Ai(p,-j/)^(p-P / W(p') (15) 
and expand the vertex in Dirac space as 

A l(k, ~P) = f^,o{k, -p) + f»,v{k, -p)l v + £ p< „ U, P u(k, ~P)lplv 
+U,5{k, -p)ib + U,u5(k, -P)lulb 

(16) 

(with an identical labeling for the decomposition in coordinate space). A lot of 
algebra |jj gives 

\(imo)[\'{im ) - p (im )] 

m B = TT- \1 — r ^ T7 \~u~- \ ( 17 ) 

P 2 !m o Ji.i '"io - ic{ 2 {tmo)A(im ) 



where 



X(im) = ^e m "°A(n) (18) 

n 

A'(im) = ^n e mno A(n) (19) 
Po (im) = -iJ2noe mno p (n) (20) 

n 

p' 2 {im) = -i^n 2 e mn V 2 (n) (21) 

n 

f 1A (im)=J2^ mX hA^y) (22) 

xy 

i 2 (im) = J> - 2y) 2 e mx f hl2 (x, y) (23) 



are all real. 



All approximate FP vertices I have seen have a complicated Pauli term with 
sizable contributions when the quark and antiquark do not sit on the same 
lattice site. As an example, Table |] shows the contribution of various fermion 
offsets to ic\ 2 (x, im) for a vertex based on the RGT we have been using in this 
paper. The normalization appears to be reasonably well saturated by fermion 
offsets over a cube. The gauge connections are very complicated. 

However, it is not clear how important the FP version of the Pauli term will 
be in spectroscopy. I therefore studied three possibilities: 
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X 


m = 0.08 


mo = 0.16 


m = 0.32 


mo = 0.64 





-0.0913 


-0.0840 


-0.0709 


-0.0527 


10 


-0.2490 


-0.2269 


-0.1883 


-0.1278 


110 


-0.2980 


-0.2710 


-0.2253 


-0.1479 


1110 


-0.1641 


-0.1493 


-0.1248 


-0.0816 


1111 


-0.0318 


-0.0291 


-0.0246 


-0.0167 


Sum in hypercube: 


-0.8343 


-0.7604 


-0.6340 


-0.4269 


Total: 


-0.8974 


-0.8182 


-0.6837 


-0.4807 



Table 2: Nonlocality of the Pauli term ic\ 2 (x,im) for actions along the RT. 

1) No Pauli term at all. This turns out to give hyperfine interactions which 
are too small. 

2) Keep only the on-site part of the Pauli term (the standard clover term) 
but choose its normalization so that tub — m . This is not a FP action. The 
gauge links will be fattened like the rest of the links in the action. 

3) Restrict the Pauli term to offsets which span a cube. For each offset, sum 
over all the minimum-length paths (with their sign factors) which contribute to 
the Pauli term. Choose the relative normalization of the terms to match the 
FP vertex and the overall normalization to fix raj = mo. Fatten the links if 
necessary. This choice (hereafter called a "full Pauli action") is a much better 
approximation to a FP action than the second choice ("clover action"), but if 
the clover action performs as well in a test, it is the action of choice. 

For m > the constraint raj = m fixes the normalization of the Pauli 
term. I find that the normalization varies roughly linearly with the bare quark 
mass. I choose (arbitrarily) to keep the same linear dependence with mo even 
for negative bare mass. This would be equivalent in the standard clover action, 
to making the size of the clover term a function of the hopping parameter k, 
rather than a function of (3. From a practical point of view the difference is 
slight: as one varies the gauge coupling in a simulation, the value of bare quark 
corresponding to a particular physical hadron mass shifts, becoming (typically) 
more negative as (3 decreases. The input coefficient of the clover term becomes 
larger as (3 decreases, so the net result is that one can as well say that the 
clover term tracks the bare quark mass, as to say that the it tracks the bare 
coupling. At zero gauge coupling FP actions do not have a mass-independent 
Pauli normalization and so the standard practice of making the normalization 
mass independent is unnatural. Presumably one could do simulations with a 
standard action, such as the Wilson-plus-clover action, tuning to equal mo, 
although it is hard to see the point of doing this as long as the dispersion relation 
is imperfect. 
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The linear parameterization reproduces mj = too to within five per cent for 
mo < 0.4. For the clover hypercubic action I could have simply set tub = mo 
for too > by inverting Eqn. [17L although I did not do that. 



3 Scaling Tests 

3.1 Survey of Actions Tested 

Most of the quenched spectroscopy has been done using a new parameteriza- 
tion of a FP gauge action for SU(3). In the Appendix I tabulate the critical 
temperature, string tension, and Sommer parameters for the gluonic action, so 
the reader can convert to his favorite scaling variable. One test has been done 
using the original parameterization of the action presented in Ref. |lq| , and 
some tests use the SU(3) gauge action of our recent work on instantons|l7|. We 
have made rough measurements of its f3 c (N t ) for deconfinement to set the scale. 

All the fermionic actions are made gauge-invariant by replacing the offsets 
by an average over the shortest distance gauge paths. For example, 

^(x)^{x + fi + i>) -> -^(x)[V ft (x)V v {x + i>) + V v (x)V fl [x + [i)]ip{x + fi + D) (24) 

where V^(x) is either one of the original links or an APE- blocked link. 

The cost of a hypercubic action per iteration step during matrix inversion 
is about 20 times as expensive as the usual Wilson action, since there are more 
neighbors and the Dirac connections are not projectors. Actions with the com- 
plicated Pauli term are about 56 times as expensive as the usual Wilson action. 
All the gauge connections are pre-computed, so there are startup and storage 
costs, as well. I used the stabilized biconjugate gradient (biCGstab) algorithm 
for matrix inversion Jig] . 

The two actions which were tested most extensively both have a hypercubic 
kinetic term and APE-blocked links with N = 7 and a — 0.3. Action A has a 
full Pauli term. It is the best approximation to a FP action I found. It used 
the very expensive gauge action of Ref. fTij ] . Action C has only the clover term 
but is otherwise identical. My scaling tests of it used the new gauge action 
presented in the Appendix. 

I tested several other actions. All the variants of tadpole improvement I 
studied had large mass renormalization. A pure hypercube action with no Pauli 
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term had a good dispersion relation at the coarsest lattice spacing, but its 
hyperfine splittings were basically identical to those of the Wilson action. 



3.2 Spectroscopy 

Lattice volumes were 4 3 x 16 at aT c = 1/2 (excessively long in the time direction, 
in retrospect), 6 3 x 16 at aT c = 1/3, and 8 3 x 16 (dangerously short) and 8 3 x 24 
(safer) at aT c = 1/4. 

The data set for Action A consists of 80 lattices at aT c = 1/2, 50 lattices 
at aT c = 1/3 and 36 8 3 x 16 lattices at aT c = 1/4. The data set for Action C 
consists of 80 lattices at aT c = 1/2 and aT c = 1/3 and 60 8 3 x 24 lattices at 
aT c = 1/4. 

The spectroscopy measurement is entirely straightforward. I gauge fixed to 
Coulomb gauge and used a Gaussian independent particle source wave function 
ip(r) = cxp(— 7r 2 ) with 7=1, 0.5, and 0.25 at aT c — 2,3,4. I used pointlikc 
sinks projected onto low momentum states. I used naive currents (ip^ip, etc.) 
for interpolating fields. The spectra appeared to be asymptotic (as shown by 
good (correlated) fits to a single exponential) beginning at t ~ 2 (at aT c = 1/2), 
3-5 (at aT c = 1/3) and 5-7 (at aT c = 1/4) and the best fits were selected using 
the old HEMCGC criterion JT§. 

My fiducial for comparison, simply because there are extensive data sets, is 
Wilson-action quenched spectroscopy. I have tried to restrict the data I used for 
comparison to lattices with the proper physical volume. I constructed my own 
aT c = 1/2 and aT c — 1/3 Wilson data sets (f3 — 5.1 and 5.54) since I could not 
find any results for these. I also ran off 40 Wilson lattices at j3 = 5.7 (N t — 4) 
to measure a dispersion relation. At that coupling my masses were within a 
standard deviation of the much superior data set of Butler, et al. ppf . 

To compare with the more standard improved actions, I also performed a 
fiducial study using the clover action (Wilson fermion action plus on-site clover 
term) in a background of Wilson action gauge fields. I tadpole-improved the 
action using u = (Trt/p/3) 1 / 4 where U p is the average plaquette. Data are at 
aT c = 1/2 and aT c = 1/3, where u Q = 0.802 and 0.844. This data set was 180 
and 60 lattices at the two couplings. 

I show first plots of m p /T c and miq/T c vs. m^jT c with Wilson fcrmions. 
These plots are scaling tests by themselves, or one can interpolate in the curves 
to fixed values of mn/Tc (equivalent to fixed quark mass) and plot the variation 
in the observable vs. aT c . Fig. || shows the rho mass and Figs. |i| shows the 
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Figure 3: m p /T c vs. m^/Tc for the Wilson action. Data are labelled with 
octagons for aT c = 1/2, diamonds for aT c — 1/3, crosses for aT c = 1/4, squares 
for aT c = 1/8, and fancy crosses for aT c = 1, 12. 



nucleon mass. Notice that the rho mass has the worst scaling violations of the 
three particles. 

We can roughly estimate the critical bare quark mass (at which the pion is 
massless) by linearly extrapolating m„ to zero in m . Fig. || shows the squared 
pion mass vs bare quark mass for action A, at f} c for N t = 2, 3, 4, and the same 
plot, but for action C, is shown in Fig. || Both actions have small bare mass 
rcnormalization. This is important from the point of view of principle because 
a true FP action would have no additive mass renormalization. It is important 
in practice because we only really know the kinetic parameters by solving the 
RG equation for positive bare mass; they must be extrapolated in some artistic 
way if one needs to go to negative bare mass. 

Another way to estimate the critical bare mass is to use or the PCAC relation 

V p • (^ 75 ?/>(0)^757^(20) = 2m 9 (^75^(0)^75^(2:)) • (25) 

If we convert to lattice operators, sum over spatial slices, and measure distance 
in the t direction, this becomes: 

Za §1 E (^75^(0)^7570^)) = 2am q Z P ^ (^(O)^^)). (26) 
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Figure 4: tun /T c vs. m v /T c for the Wilson action. Data are labelled as in Fig. 
11. 



I follow |21J by fitting the pseudoscalar source-pseudoscalar sink to 

P(t) = Z(exp(-m 7r t) + exp(-m Tr (7V t - t))) 
and the pseudoscalar source-axial sink to 
Z Qui 

Mt) = if q -Z(exp{-( m7I t) - exp(-m T (N t - t))). 



(27) 



(28) 



to extract m q . There are many other possibilities for defining an axial current 
and for defining the derivative operator. I only use the naive (pointlike) currents. 
I do not know the Z— factors, but for finding the value of mg that does not 
matter. Extrapolating or m q linearly in mo ignores all the well-known 
problems associated with extracting quark masses from lattice data [p2[ , but 
the procedure is perfectly adequate to distinguish a small quark mass from a 
large one. The quark masses are shown in Figs. | and |. 



We find for action A that 



-0.42, -0.20, and -0.14 at aT c = 1/2, 



1/3, and 1/4, respectively. For action C, the corresponding numbers are mg = 
—0.36, —0.16, and —0.095. These numbers should be compared to the analogous 
quantities for Wilson fermions, using mg = l/(2/c c ) — 4: -1.58 at f3 = 5.1 
(aT c = 1/2), -1.04 at (3 = 5.7 (aT c = 4), and still -0.70 at j3 = 6.3 @. In the 
latter case tadpole improved perturbation theory can explain most of the mass 
shift. 



Now for scaling tests. I compare m p /T c and mjy/T c vs. m^/T c for actions 
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Figure 5: Bare squared pion mass (upper set of curves) and quark mass from 
Eq. 40 (lower set of curves) vs bare quark mass for action A, at (3 C for N t = 2 
(squares), 3 (diamonds), 4 (crosses). 



-0.4 



-0.2 







0.0 



0.2 



Figure 6: Bare squared pion mass (upper set of curves) and quark mass from 
Eq. 40 (lower set of curves) vs bare quark mass for action C, at (3 C for N t = 2 
(squares), 3 (diamonds), 4 (crosses). 
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Figure 7: m p /T c vs. m^/T c for action A. 



A and C in Figs. @, §, g, and 0. 

We can compare scaling violations in hyperfme splittings by interpolating 
our data to fixed Tr/p mass ratios and plotting the N/p mass ratio vs. m p a. I 
do this at four ir/p mass ratios, 0.80 and 0.70, in Fig. O. In these figure the 
diamonds are Wilson action data in lattices of fixed physical size (4 3 at (3 = 5.1, 
6 3 at f3 = 5.54, 8 3 at (3 = 5.7 |o), 16 3 at /3 = 6.0 |§ 24 3 at /3 = 6.3 @) and 
the crosses are data in various larger lattices: 16 3 and 24 3 at (3 — 5.7 and 32 3 at 
13 = 6.17 @, 24 3 at (3 = 6.0 (23). When they are present the data points from 
larger lattices illustrate the danger of performing scaling tests with data from 
different volumes. The bursts are from the nonperturbatively improved clover 
action of Ref. |^5| and the fancy diamonds are the TI clover action. The other 
plotting symbols show our test actions A and C. 

We give a tables of masses from Action C in Tables 0, ||, and ||, and for 
Action A in Tables |, 0, and |. 



3.3 Dispersion Relations 



There are two ways to look at dispersion relations. The simplest is to plot E(p) 
the energy of the state produced with spatial momentum p, as a function of \p\ . 
The result for action C at bare mass 0.15 is compared to the free dispersion 
relation at aT c = 1/2 in Fig. [12[ All of the test actions I have studied have 
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Figure 8: mjv /T c vs. m^/Tc for action A. 
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Figure 9: m p /T c vs. ra„/T c for action C. 
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Figure 10: itin/T c vs. m n /T c for action C. 



good dispersion relations even at aT c = 1/2. I believe that is a generic feature 
of the hypercubic kinetic term. 

The signal to noise ratio for the nonzero momentum meson channels dies 
away at large t like exp(— (E(p) — m w )t). This means that large statistics are 
required to go to small quark mass or to high p. However, it is possible to extract 
a fitted c 2 e ^ = (E(p) 2 — rn 2 )/p 2 for the lowest nonzero momentum mode, for 
larger masses. This was done by performing a 4-parameter correlated fit to a 
pair of single exponentials, one for the p — (0, 0, 0) mode and the other the 
p = (1,0,0) mode. I compare my results from action C and from the Wilson 
and clover actions in Figs. |l3| and |lj. Hadron masses are again scaled by T c to 
allow the display of several lattice spacings at once. Action C has ~ 1 for 
all observed hadrons even at p = ttT c (860 MeV/c) at aT c — 1/2. 



3.4 Summary 

It appears that these actions are members of a family of actions which show 
improved scaling, even at j3 c {N t = 2), about 0.36 fm lattice spacing. 

The hypercubic actions have much better dispersion relations than either 
the clover or Wilson action. They share this improvement with the D234 family 
of actions [^6| and with the Hamber-Wu [^?J action as tested in Ref . |^8| . 
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Figure 11: A scaling test for new actions: (squares action C, octagons action 
A) vs. Wilson actions on lattices of fixed physical size (diamonds) and larger 
volumes (crosses), and the nonperturbatively improved (bursts) and tadpole 
improved (fancy crosses) clover actions. Data are interpolated to n/p = 0.84 
(a), 0.80 (b), 0.75 (c) and 0.70 (d). 
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arriq 


PS 


V 


N 


A 


0.30 


1.883( 8) 


2.066( 8) 


3.332(17) 


3.450(21) 


0.15 


1.623( 9) 


1.861(10) 


3.002(20) 


3.156(25) 


0.10 


1.535( 9) 


1.796(11) 


2.891(22) 


3.061(27) 


0.05 


1.446( 9) 


1.733(12) 


2.781(24) 


2.967(30) 


0.00 


1.363( 9) 


1.672(14) 


2.670(27) 


2.875(33) 


-0.05 


1.251(10) 


1.575(20) 


2.526(37) 


2.770(46) 


-0.10 


1.159(10) 


1.551(20) 


2.432(39) 


2.687(43) 


-0.15 


1.039(12) 


1.491(26) 


2.303(54) 


2.591(49) 


-0.20 


0.899(16) 


1.429(35) 


2.147(70) 


2.505(64) 


Tabic 3: Table of best-fit 


masses, action C, (3 = 3 


am q 


PS 


V 


N 


A 


0.15 


1.156( 7) 


1.363(11) 


2.077(21) 


2.258(28) 


0.10 


1.066(15) 


1.328(20) 


2.090(33) 


2.277(51) 


0.05 


0.935( 8) 


1.218(15) 


1.816(30) 


2.049(31) 


0.00 


0.825( 9) 


1.153(19) 


1.688(38) 


1.959(41) 


-0.05 


0.689(10) 


1.124(18) 


1.603(27) 


1.941(33) 



.092 (aT c 



Table 4: Table of best-fit masses, action C, [3 — 3.50 (aT c 



am q 


PS 


V 


N 


A 


0.15 


0.943( 4) 


1.118( 9) 


1.722(10) 


1.853(14) 


0.10 


0.836( 4) 


1.044(11) 


1.592(12) 


1.742(17) 


0.07 


0.768( 5) 


0.997(12) 


1.514(14) 


1.678(19) 


0.03 


0.672( 5) 


0.934(15) 


1.410(18) 


1.598(22) 


Table 5: Table of best-fit masses, action C, (i = 


am q 


PS 


V 


N 


A 


0.15 


1.816(11) 


1.972(18) 


3.274(36) 


3.355(44) 


0.10 


1.727(11) 


1.903(19) 


3.161(37) 


3.309(49) 


0.05 


1.637(11) 


1.835(21) 


3.047(39) 


3.237(47) 


0.00 


1.546(11) 


1.769(22) 


2.930(43) 


3.156(51) 


-0.10 


1.358(10) 


1.582(20) 


2.596(44) 


2.825(43) 


-0.15 


1.246(10) 


1.502(24) 


2.452(44) 


2.718(49) 


-0.20 


1.123(11) 


1.417(29) 


2.297(45) 


2.605(57) 


-0.25 


0.982(13) 


1.322(37) 


2.155(42) 


2.483(71) 



3.70 (aT c 



Table 6: Table of best-fit masses, action A, (3 = 2.38 (aT c 
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Figure 12: Dispersion relation for heavy hadrons at aT c = 1/2 (a ~ 0.36 fm) 
from action C. The curves are the continuum dispersion relation for the appro- 
priate (measured) hadron mass. 



am q 


PS 


V 


N 


A 


0.10 
0.05 
0.00 
-0.05 


1.179( 8) 
1.072( 8) 
0.956( 9) 
0.829(10) 


1.353(11) 
1.275(12) 
1.194(13) 
1.113(16) 


2.145(19) 
2.010(21) 
1.870(24) 
1.720(28) 


2.245(24) 
2.115(27) 
1.989(32) 
1.854(40) 



Table 7: Table of best-fit masses, action A, (3 = 2.85 {aT c = 1/3). 



However, the hypercubic actions tested here seem to produce about the same 
level of improvement in hyperfine splittings as the clover action, at heavier quark 
masses. Leaving out the Pauli term gives noticeable scaling violations with a too- 
large N/p ratio; probably one needs to keep some kind of explicit Pauli/clover 
term in the lattice action to boost the hyperfine splittings. 

The hyperfine splittings show worse scaling violations than the dispersion 
relation. Controlling and approximating the quark anomalous magnetic moment 
is the most difficult part of the construction of a FP action, and that may be 
the source of the scale violations. 

Of course, there is still the possibility that all the actions tested here are 
missing some other common physics ingredient, which is responsible for scaling 
of the hyperfine splittings. 

Since the p meson is the particle which shows the largest scaling violations, 
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Figure 13: Squared speed of light vs. liadron mass in units of T c , for (a) 
pseudoscalars, (b) vectors) and (c) protons, from action C. Octagons, crosses, 
and diamonds label aT c = 1/2, 1/3 and 1/4. 



am q 


PS 


V 


N 


A 


0.15 
0.10 
0.05 
0.00 


1.100( 5) 
0.993( 6) 
0.880( 6) 
0.758( 7) 


1.228( 8) 
1.146( 9) 
1.066(11) 
1.000(12) 


1.903(20) 
1.783(18) 
1.640(20) 
1.504(25) 


1.998(27) 
1.875(30) 
1.788(26) 
1.662(29) 



Table 8: Table of best-fit masses, action A, (3 = 3.05 {aT c = 1/4). 
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Figure 14: Squared speed of light vs hadron mass in units of T c , for (a) 
pseudoscalars, (b) vectors) and (c) protons, from the Wilson and clover actions. 
Octagons, crosses, and diamonds label aT c = 1/2, 1/3 and 1/4 for the Wilson 
action, and for the clover action the labels are plusses for aT c = 1/4 and fancy 
crosses for aT c = 1/3. 
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the best way to quantify improvement is by taking "sections" of the m p /T c and 
toat/Tc vs. m^/Tc plots and displaying them vs. aT c at fixed m w /T c in Fig. 
p"5| . The improved actions at aT c = 1/2 seem to show the same level of scaling 
violations as Wilson data at a lattice spacing a factor of 3 smaller. The smaller 
lattice spacing data seem to pick up about a factor of two improvement in lattice 
spacing, although the uncertainty in the data is larger. Quenched calculations 
are thought to scale in difficulty like 1/a 6 ; the cost of action C is about a factor 
of 20 compared to the Wilson action. 

My data by themselves do not suggest a unique way to extrapolate to a = 
0. FP actions are classically perfect with no a™ scaling violations for any n. 
Approximate FP actions generally have discretization errors at all orders in a, 
though the coefficients of any order in a are typically much smaller than the 
corresponding coefficient in an improved action (for n > N). 



4 Conclusions 



The important ingredients of these actions which contribute to their improved 
scaling behavior are the hypercubic kinetic term and the lattice anomalous mag- 
netic moment term. The very small additive renormalization of the bare lattice 
mass is due to the use of fermion- gauge field couplings which are insensitive to 
the short-distance fluctuations of the gauge fields. 

The specific implementation of these ideas in the actions I have tested in- 
volves many arbitrary choices. 1 believe that essentially every choice I made 
for a particular parameterization could be replaced by another choice, which 
would give an action which would have the same quality of scaling violations. 
Some changes should add some additional good feature. For example, it might 
be possible to find a parameterization of a fat link which would lend itself to 
simple perturbation theory calculations. 

There are several obvious extensions of this work. The first involves the 
kinetic term: It would be useful to find a parameterization of the kinetic term 
which extends from zero mass to very large mass. It would also be interesting 
to develop a fat link parameterization of the action which could be efficiently 
incorporated into one of the standard algorithms for dynamical fermions. Next, 
is there a better parameterization of the anomalous magnetic moment term 
which might improve scaling? 

Fixed point actions have many desirable formal properties JToj ] : they include 
scale invariant instanton solutions, the index theorem, an absence of exceptional 
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Figure 15: Variation of m p /T c and /T c vs. aT c at fixed m^/T c = 2.75 ((a) 
and (b)) and 3.15 ((c) and (d)), for actions A (octagons), C (squares), Wilson 
action (diamonds), and TI clover action (fancy crosses). 
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configurations, and a remnant of chiral symmetry |2^, [30|. These properties may 
not be present in an action which is a bad approximation to a FP action. I do 
not know how well they are satisfied by these actions (other than the apparent 
small renormalization of the quark mass.) 

Indeed, the particular choice of a free field action which I made was motivated 
only by the locality and spectral properties of the free action. No attempt was 
made to optimize the chiral properties of the approximate action. This is clearly 
the outstanding problem for future study. 



A Pure Gauge Actions 



This work used a new few-parameter approximation to an FP action for SU(3) 
gauge theory. Using it, isolated instanton configurations have constant actions 
to within 1.5 per cent. Like the action of Ref. |16| ], it is a superposition of 
powers of the plaquette and the perimeter-6 "twisted" link (x, y, z, —x, —y, —z). 
Like the action of Ref. jlTj it includes a constant term. It is designed to be 
used for couplings such that the lattice spacing is aT c ~ 1/3 or 1/4 to 1/8 or 
so. Explicitly 

S(V) = Co + ±- 5>i(C)(tfc - Tr(Vc)) + c 2 (C)(N c - Tr(V c )) 2 + ... (A.l) 

c c 

with coefficients tabulated in Table pi 



Table 9: Couplings of the few-parameter approximate FP action. 



operator 


Cl 


C2 


C3 


C4 


c = -2.517 


Cplaq 


3.248 


-1.580 


.1257 


.0576 




C^—link 


-.2810 


.0051 


.00490 


-.0096 





This action costs about a factor of 7 times the usual Wilson plaquette action 
to simulate. 

I have measured the critical couplings for the deconfinement transition at 
aT c = 1/2, 1/3, and 1/4. The critical couplings on the measured spatial volumes 
and my extrapolation to infinite volume are shown in Table [l^. I have also 
measured the string tension from Wilson loops at these values of the coupling, 
on 8 4 lattices (at (3 = 3.092) and 12 4 lattices for the other couplings. The data 
was fitted to a static potential V(r) of the form 

V(r) = V + ar - E/r (A.2) 



26 



using the techniques of Ref. The fit to the largest lattice spacing data 

is very difficult. The signal from large r is not good, and there is very little 
left of the Coulomb part of the potential due to the coarseness of the lattice. 
Nevertheless, I present the string tension and the Sommer p2j parameter r 
{r^dV {tq) / dr — —1.65) in Table |l^. We see scaling within errors for both these 
parameters (vs. T c ) at aT c = 1/3 and 1/4. There is a ten per cent scaling 
violation at aT c = 1/2. The asymptotic value inferred from large scale Wilson 
simulations §| is s/a/T c = 1.600(11). 

Finally, I show a plot of V(r)/T c vs rT c for the three lattice spacings in Fig. 
[L6[ The overall vertical shift in the potentials is not physical, but it allows the 
reader to separate the different data sets by eye. 



Table 10: 


Critical couplings at finite volume and extrapolated to infinite volume 


for the FP action with parameters 


in Tabic 1. 




volume 


N t = 2 


N t =3 


N t = A 


4 3 


3.025(25) 






6 3 


3.06(1) 


3.47(1) 




8 3 


3.08(1) 


3.49(1) 


3.67(1) 


10 3 


3.085(5) 


3.50(1) 


3.69(1) 


12 3 






3.69(1) 


infinite 


3.092(7) 


3.50(1) 


3.70(1) 


T c /A 


8.67 


8.96 


8.33 


•> 

a a 


0.56(5) 


.302(16) 


.164(3) 




1.50(7) 


1.65(2) 


1.62(2) 


r /a 


1.3(2) 


2.19(2) 


2.93(1) 


r T c 


.65(10) 


.730(7) 


.733(3) 
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